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Rare nonleptonic f2~ — > H°(1530)7r~ decay branching ratio is estimated by means of the QCD 
enhanced effective weak Hamiltonian supplemented by the SU(3) Skyrme model used to estimate 
the nonperturbative matrix elements. Using mean values for experimental input parameters and 
the Skyrme charge e = 4.75 we obtain the rate which is in a good agreement with data. 
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It is well known that the nonleptonic weak decays of 
baryons can be reasonably well described in the frame- 
work of the Standard Model 0, Recently both s-, 
p-wave nonleptonic hyperon and p-wave £l~ decay am- 
plitudes were quite successfully reproduced by the SU(3) 
extended Skyrme model with the QCD enhanced effec- 
tive weak Hamiltonian 0, 0| ■ The decay amplitudes were 
described through the current-algebra commutator, the 
ground-state baryon pole terms and factorizable contri- 
butions. The nonperturbative quantities, i.e. the baryon 
six dimensional operator matrix elements, were deter- 
mined using the SU(3) extended Skyrme model. For the 
s-wavc hyperon decay amplitudes correct relative signs 
and absolute magnitudes were obtained. For the p- waves 
all relative signs were correct, with their relative mag- 
nitudes roughly following the experimental data. One 
is thus faced with an obvious question: could an analo- 
gous approach work equally well for the rare nonleptonic 
il~ — > Sgy 2+ (1530)7T~ H weak decay? Such a question 
should be considered in connection with the measure- 
ments of the fl~ lifetime and branching ratios |(|: 
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In this report our goal is to test whether the effective 
weak Hamiltonian and the SU(3) Skyrme model are able 
to predict the rare nonleptonic f2~ — ► Z* a ir~ decay fol- 
lowing the method in 0, [j] . To this end we shall employ 
the Standard Model effective Hamiltonian and the min- 
imal number of couplings concept of the Skyrme model 
to estimate the nonperturbative matrix elements of the 
4-quark operators 0, 0, f° r decuplet-decuplet transi- 
tions. Throughout this report we use the arctan ansatz 
for the Skyrme profile function F(r) Q, which allows to 
calculate the pertinent overlap integrals analytically with 
accuracy of the order of ~1 % with respect to the exact 
numerical results. 

The starting point in such an analysis is the effective 
weak Hamiltonian in the form of the current <X> current 



interaction, enhanced by QCD, 

Hgf(AS - 1) = V2G F V; d V us <*Oi, (3) 

where Gf is the Fermi constant and V* d V us are the 
Cabbibo-Kobayashi-Maskawa matrix elements. The Oi 
are the 4-quark operators and the c, factors are the 
QCD-short distance Wilson coefficients 0,0: ci = 
-1.90-0.61C, c 2 = 0.14 + 0.020C, c 3 = c 4 /5, c 4 = 
0.49+0.005C, with C = V* d V ts /V* d V us . For the purpose of 
this work we neglect the so called Penguin operators since 
their contributions are proven to be small Q. Without 
QCD corrections, the Wilson coefficients would have the 
following values: c\ = —1, c 2 = 1/5, C3 = 2/15, C4 = 
2/3. In this paper we consider both possibilities and 
compare the resulting rates. 

The techniques used to describe nonleptonic Q ~ decays 
(in this work we have only the 3/2+ — * 3/2+ + 0~ reac- 
tion) are known as a modified current-algebra (CA) ap- 
proach. The general form of the decay amplitude reads: 



{^{q)_B'{p')\H^f\B{p)) = 



(4) 



The W(p) denotes the Rarita-Schwinger spinor. The 
parity-violating amplitudes A correspond to the s-wave 
and parity-conserving amplitudes B correspond to the p- 
wave f2 _ decays, respectively. Since the decay of S*°7r~ 
is strongly limited by phase space (momentum transfer ~ 
1 McV at the peak value of the S* mass), we will neglect 
the amplitudes C and T> (d- and f- waves). The decay 
probability T(tt- /2+ (p) -> E*° /2+ (p') + tt". (q)) is: 
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Here ft is at rest, m/ denotes the final baryon mass and 
rrirf, is the mass of the emitted meson. 
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We calculate the decay amplitudes by using the so- 
called tree-diagram approximation at the particle level, 
i.e. factorizable and pole diagrams plus the commutator 
term. The amplitudes Q receive contributions from the 
commutator term and the pole diagram, respectively: 
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However, another state rather strongly coupled to S*7r 
and with a mass close to the il~ exists, the 5(1820, J p = 
3/2 ~) resonance 0, whose mass fits nicely to the Gell- 
Man-Okubo formula for an octet of 3/2~ baryons. There- 
fore, there is a pole term in the s-wave amplitudes, A-p: 



A, 



win - ms* 
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The commutator and the baryon-pole amplitudes contain 
weak matrix elements defined as 



■^ G ^v us ^ ]ci0rm ), (11) 



where the important parts are the Wilson coefficients and 
the 4-quark operator matrix elements. 

The factorizable contributions to s- and p-waves, are 
calculated by inserting the vacuum states; it is therefore a 
factorized product of two current matrix elements, where 
the decuplet-decuplct matrix element of the vector and 
the axial-vector currents reads: 



<H*(p')| 



|n _ (p)) = sg^(p) 7 ^Jw>')(i2) 



Summing over all factorizable contributions gives the fol- 
lowing expressions for the amplitudes: 
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V: d V us (maTms*)U (13) 



x g% n [ci - 2(c 2 + c 3 + c 4 )] 



SU(3) extended Skyrme model action is C — C a + Csk + 
Csb + £wz, where C a , Csk, £.sb, and Cwz denote the 
cr-model, Skyrme, symmetry breakin g (SB) and Wess- 
Zumino (WZ) terms, respectively ^ . 1 171 Hal . Extension 
of the model to the strange sector 13L Il5l fl6| is done 
by an isospin embedding of the static hedgehog ansatz 
into an SU(3) matrix, which is a subject of a time de- 
pendent rotation U(r,t) = A(t)U (r) (t) by a collective 
coordinate matrix A(t) G SU(3). The generalized veloc- 
ities are defined by A^(t)A(t) = and the 
profile function is interpreted as a chiral angle that pa- 
rameterizes the soliton. The collective coordinates a a are 
canonically quantized to generate the states that possess 
the quantum numbers of the physical strange baryons. 
In order to account for a non-zero strange quark mass 
the appropriate symmetry breaking terms should be in- 
cluded. 

In this work we will use the SU(3) extended action 
C with the following set of parameters [l^, x = 36.97, 
P' = -28.6 MeV 2 , 6' = 4.12 x 10 7 MeV 4 , determined 
from the masses and decay constants of the pseudo-scalar 
mesons. The x term is responsible for the baryon mass 
splittings and the admixture of higher representations in 
the baryon wave functions. However, we use the SU(3) 
symmetric baryon wave functions in the spirit of the per- 
turbative approach to SB. Indeed, we have shown in 
0, 0, that the WZ and SB contributions through the 
weak operator matrix elements are small and introduce 
few % uncertainty of the decay amplitudes dominated by 
the Skyrme term which scales like 1/e. 

Since the coupling f n is equal to its experimental value, 
for the evaluation of rare nonleptonic decay rate the 
only remaining free parameter, as in 0,0, is the Skyrme 
charge e. It has been shown in Figure 12 of 0] that for 
4.0 < e < 5.0 the mass spectrums of 8, 10 and 10-plcts 
are reasonably well described. This is the reason we are 
using that particular range of e further on. 

The 4-quark operator matrix element contribution to 
the commutator |JSJ and to the baryon pole term 0, 
evaluated as a function of the Skyrme charge e 0, in 
units of [10~ 8 GeV], is 



where the g^iy) represents the form-factor of the spatial 
component of the axial- vector (vector) current. 
The total theoretical amplitudes are: 

Ath(ml) - A Comm (0) + A v (ml) + A s (ml), (14) 
B th (ml) = Bvim^ + Bsiml), (15) 

were the relative signs between commutator, pole and 
factorizable contributions are determined via SU(3) and 
the generalized Goldberger-Treiman relation. 

In order to estimate the 4-quark operator matrix el- 
ements entering (fill) , we take the Skyrme model where 
baryons emerge as soliton co nfig urations of the field U 
of pseudo-scalar mesons [H Hi E Q III The 
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QCD off 


QCD on 


e 


-1.86 


-3.52 


4.00 


-1.31 


-2.75 


4.75 


-1.22 


-2.55 


5.00 



(16) 



The size of the strong coupling from pole term 10, 
determined via generalized Goldberger-Treiman relation, 



= — ^—r-- — ^^9a 
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9-11 g p A n , (17) 



is close to the estimate given in 

The strong coupling gs-^s'^, which enters in the cal- 
culation of A-p, we extract from experimental value for 
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partial width — > H*tt). Using strong effective La- 
grangian for spin 3/2~ baryon 5** 

£ int (S** -> S*tt) = g 3 **3*.W„(p')W»(p)Mq), (18) 

where VV M denotes the 3/2~ Rarita-Schwinger spinor, we 
obtain for S** at rest, 



9i 



■(S' + ma.) (19) 

3b7rm=:»mH** 
x [2E a + 2E'm s , + 5ml,] . 

The E' and |q| are given in JfjJ) and J7J. From experiment 
we have r tot (S**) = (24 ± 6) MeV and 

r(s** -> s*7r)/r tot (s**) = (0.30 ± 0.15). (20) 

The generalized Goldberger-Treiman relation applied 
to ijl"U|) gives the product &H**n<?E**E* which docs not 
depend on the phase of the exchanged S** in the pole 
diagram. It has been found in [2jj that the sign of an 
axial- vector matrix clement, between particles of the op- 
posite parity, is negative. Than from l|19|) and (|20fl we 
find <7h«h*7t = — (0.48 ± 0.14), a small negative value 
with large uncertainty. 

The weak matrix elements of the same spin but oppo- 
site parity baryons are found to be about the same as the 
ones including ground states only [2(|, i.e. &h»*si = a~*Q- 
This together with © and (pll gives 



A v = -(0.44 ±0.12) A 



Comjn i 



(21) 



up to the sign, similar to the conclusion in [19|. Since 
the resonance S** is known with quite a large error, the 
same error is also present in <|21ll . 

We proceed with the computation of the vector and 
axial- vector current form-factors ffy(5) > i n terms of gyn^ , 
via SU(3) Clcbsh-Gordan coefficients. Using from exper- 
iment, q'fVqF" = 1.26, we find an agreement with the 
SU(6) result [lj, as we should: 
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Factorizablc amplitudes A$ and B$ arc very slowly 
changing functions, via the q^ n / 'g v n , with respect to the 
Skyrme charge e, (for details see Fig's 2-4 in (ljjV For 
QCD off/on we found from <P|l and Q : 



Af 
Bf 



-5.65 x 10 
-90.1 x 10 



-8 a on 
-8 non 

> B s 



-6.32 x 10" 8 , 
-100.8 x 10~ 8 . (23) 



The mean values of theoretical amplitudes Acomm, A-p 
and B-p for QCD off /on, as a functions of the Skyrme 
charge e, are presented in Table U 

Rare nonlcptonic two-body decay partial width for- 
mulae 

r th (fT -> E*°tt-) = {9.000617|A 5 + Acomm + A v \ 2 
+ 0.000148|B,s + B V \ 2 } 2.697 x 10" 4 GcV, (24) 



TABLE I: The amplitudes A Comm , A v , B v in units [10 8 ] 
contributing to the rare nonleptonic two-body Q.~ decay rate. 





QCD off 


QCD on 


e 


Acomm Av B-p 


Acomm A-p Bp 


4.00 
4.75 
5.00 


14.14 -6.22 155.6 
9.96 -4.38 109.6 
9.28 -4.08 102.1 


26.76 -11.77 294.5 
20.91 -9.20 230.1 
19.39 -8.53 213.4 



and eqs. $ to Ip^. in units of [10~ 18 GeV], gives 

QCD off QCD on 



+2.65 
1.09 
+0.38 
0.01 

06+ 54 



1.19 
0.01 



17.81 



+ 16.57 
10.44 

6.84± 8 f 7 
4.83±|;gg 



4.00 
4.75 
5.00 



(25) 



which, compared with the experiment 



r exp (rr 2*%-) = 5.1 



-4.1 
-1.6 



x 10" 18 GeV, (26) 



shows the following: 

(a) In this dynamical scheme framework, contrary to 
the nonleptonic hyperon and £1~ decays 0,0 > the fac- 
torizable contributions turn out to be very important for 
the rare fi~ — ► S*°7r~ nonleptonic decay. The opposite 
signs between the commutator and pole term (|21fl and 
the factorizable contributions (|23|l becomes an essential 
feature, of our dynamical scheme, leading to the internal 
cancellation within the A and B amplitudes and bringing 
theoretical estimate l|25f) closer to the experiment. 

(b) From theoretical rate 1|25() it is clear that in this 
dynamical scheme the QCD enhancement is crucial. 

(c) Important feature of the particular rare decay mode 
— > S*°7r - lies in the fact that it is strongly limited 

by phase space and that the decay transition occurs at 
threshold. Thus, 7r~ in the final state is almost at rest. 

(d) First consequence of (c) is the emitted pion is re- 
ally soft and the so called soft pion limit theorem is very 
well satisfied. This is the reason for the commutator 
term dominance in the s-wave amplitude A. Second con- 
sequence of (c) is an overall phase space enhancement, 
by factor of 6.1 x 10 4 , for the s- waves. Owing to this, 
altho the p-wave amplitudes are more than one order of 
magnitude larger than the corresponding s-waves, the to- 
tal contribution of B amplitudes to the rate (|25|l is very 
small. 

(e) Inspection of the result l|25() shows the importance 
of the e dependence and that e = 4.75, with QCD cor- 
rections switched on, gives the central fit for the rare 
n~ — > S*°7r~ nonleptonic decay. Applying the same e 
and QCD switched on to the £l~ nonlcptonic decays from 
0, we first find a AH o = 4.85 x 10" 8 GeV. This together 
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with the matrix element Gte*-n- = —2.75 x 10 -8 GeV 
and factorizable contributions from Table 1 in pro- 
duces the fl~ nonlcptonic decay amplitudes: Bthi^x) = 
6.02, B th (SlZ) = 1-48 and B th (^o) = 1-29, in units of 
[10~ 6 GeV -1 ], which are closer to the conclusions in [2l| 
and in better agreement with experiment than the ones 
obtained in The same should hold for the p-wave 
nonleptonic hyperon decays Q. 

Altho our dynamical scheme is the same as the one in 
|l9j |. situation has changed during the past twenty years. 
The important improvements happened in the areas of 
QCD corrections to the effective weak Hamiltonian and 
in the further development of the Skyrme model. The 
latter becomes a good candidate for the correct descrip- 
tion of the higher SU(3) representation multiplets [lUE^ 
and a good tool for evaluation of the nonperturbative 
quantities like the 4-quark operator matrix elements be- 
tween different baryon states. Comparing the result 1|25|) 
with Table 1 from [J^, we conclude that the dynamical 
scheme is working well in both cases. In this work we 
are using well accepted result for the NLO computations 
of the QCD corrections According to Q the errors 
generated by different choice of rcnormalization scheme 
and scale are less than 10%. Overall QCD corrections are 
more stabile and smaller today; penguins are about one 
order of magnitude smaller and also c\ ~ —3 — > —1.9 0. 
Skyrme model provide us with 4-quark operator matrix 
elements, in different way than the MIT bag model did 
in P, Il9| , producing rough agreement with measurement 
(|26J) . Finally, according to (e), as a bonus we obtain bet- 
ter agreement with experiment for the nonlcptonic f2~ 
and the p-wave hyperon decays. 

Obviously, not all details are under full control. For ex- 
ample m s corrections are neglected and the uncertainties 
contained in the input parameters (j20(l arc unfortunately 
large. Nevertheless, the QCD-corrected weak Hamilto- 
nian H^ff, supplemented by the Skyrme model, for the 
mean values of the input parameters, lead to the correct 
description of the rare fl~ — ► S°(1530)7r _ decay. 
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